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In a recent paper, Liddle and Urena-Lopez suggested that to have a unified model of inflation 
and dark matter is imperative to have a proper reheating process where part of the inflaton field 
remains. In this paper I propose a model where this is possible. I found that, incorporating the effect 
of plasma masses generated by the inflaton products, enable us to stop the process. A numerical 
estimated model is presented. 

PACS numbers: 98.80.Cq, 95.35.+d 



o 
o 

(N 
u 

< 

OO 



(N 
> 

(N 

o 
o 

Or 

6 



I. INTRODUCTION 

In a recent paper [l[ , Liddle and Urena-Lopez studied 
the conditions under which we can have a unified de- 
scription of inflation Q , dark energy [!| and dark matter 
0| . The key ingredient required is to have a reheating 
process where not all the inflaton energy density decays 
into radiation. The crucial role of reheating in these 
type of models was found earlier in 5], where models 
of quintessential inflation were proposed. The authors 
of [1| found that the standard reheating mechanism [2] 
can not be used for this task, because in this scenario 
the scalar field decays completely. They also studied the 
preheating scenario where it is known that the field 
decay can be incomplete, however as was discussed in [l[ 
the simplest model based on a theory with a quadratic 
scalar field turns out inconsistent with observations. 

In this work I revisited this issue. Specifically I am 
interested in a unified description of inflation and dark 
matter. The key problem is to consider an appropriate 
phase of reheating that both serve as a bridge between 
inflation and the standard model of cosmology, and also 
give us an observationally consistent amount of dark mat- 
ter. I found that even using the standard (perturbative) 
reheating mecanism, we can obtain a partially decaying 
inflaton field. However, this observation does not help to 
build the model by itself, because for a large initial am- 
plitude of the inflaton field, which is ~ M p at the end of 
inflation, we can not neglect the parametric resonance ef- 
fects [1] . So, we have to consider the preheating scenario 
, in which the first stages of reheating occur in a regime 
of a broad parametric resonance, then the resonance be- 
comes narrow, and the last stage can be described by the 
elementary theory of reheating However, this time 
the last stage is not so simple as the standard one, be- 
cause the particles created at the first stages affects the 
evolution of the inflaton field. One way to consider this, 
is take into account the generation of plasma masses for 
the inflaton decays products [|| , which can stop the par- 
ticle creation process, doing it kinematically forbidden. 
In section II we describe the basics of the model we stud- 



ied. The observational constraint is discussed in section 
III, and applied to the model in section IV where is made 
manifest the problem. The solution is described in sec- 
tion V where the plasma masses are taking into account. 



II. THE MODEL 

Let us assume a theory with a quadratic scalar poten- 
tial V{4>) — Vq + m 2 <f) 2 /2. Here Vq is a small positive 
constant needed to explain dark energy Q . During infla- 
tion, the friction produced by the expansion of the uni- 
verse makes the field evolve slowly towards its vacuum, 
e.g. m <C H. In this case the equations controlling the 
evolution are 



H 
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3H(f> + m 2 



0. 



(1) 



Inflation last until the kinetic energy of the field equals 
the potential energy cf 2 ~ V(<f>). As is well known 



the field at the end of inflation takes the value 



Of course the end of inflation coincides with the condi- 
tion m ~ H, as can be seen clearly from the Friedman 
equation. After inflation the universe enters into the re- 
heating phase, the process where almost all the particles 
in the universe were created. During this phase the scalar 
field continues rolling down the hill of the potential to- 
wards its minimum and starts to oscillate around it. In 
numerical estimates one realizes that during the first os- 
cillations the expansion of the universe is still important 
0] and the amplitude of the field falls down very quickly 
(see Figure 1). After the first oscillation, the amplitude 
reaches the value 0.04M P , indicating that the expansion 
is still important. Later the scalar field enters into the 
oscillatory regime where the amplitude decreases slowly 
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sin(mt) = sin(mi). 



(2) 
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It is during this stage where the average over many os- 
cillations of the scalar field can be described as non- 
relativistic matter ~ a~ 3 . Reheating occurs when 
the amplitude of the field decreases more rapidly than 
Eq.©. Historically this process was studied first intro- 
ducing an ad-hoc term in the equation of motion for the 
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FIG. 1: Oscillations of the inflaton field after inflation. The 
value of the scalar field is measured in units of M p , and time 
in units of mT 1 . 



field 



(j) + 3H<j> + T<t> + m 2 <f> = 0, 



(3) 



where T is the rate of particle decay of the scalar field into 
other particles. For example, if this scalar field decays in 
two scalar fields <fi — > \Xi tne ra te is 



5 V 2 



(4) 



where I am assuming a coupling g 2 ax 2 <t>- The decay 
products of the inflaton field are ultrarelativistic (m 3> 
m x ), and their energy density decreases due to the ex- 
pansion of the universe much faster than the energy of 
the oscillating field (f>. In this case reheating may ends 
when the Hubble parameter H ~ 2/3t becomes smaller 
than r. However, for a <C where $ is the amplitude 
of the oscillations, we can write 0]: T ~ <7 4 $ 2 /87tto, then 
because $ 2 decreases as i -2 (see Eq.@) in the expanding 
universe, whereas the Hubble parameter decreases only 
as t , the decay rate never catches up with the expan- 
sion of the universe, and reheating never completes. 

So even in the perturbative regime of reheating, we 
have a mechanism to stop the process leaving part of the 
oscillating field (f> decoupled and behaving as dark matter. 
However, as we can see in the next sections, this model 
by itself can not work, because it gives a result still in 
contradiction with the observational constraints. 



III. OBSERVATIONAL CONSTRAINTS 

In this section we derive a constraint for the initial 
amplitude of the scalar field <f> oscillations, which can be 



identified as the dark matter component of the model. 
The idea is to compare the theoretical results with an 
observational constraint. 

As we saw in the previous section, during inflation the 
scalar mass satisfy m <^ H , meanwhile for the quadratic 
potential the condition m ~ H marks the end of inflation. 
To recover the standard dark matter scenario, the scalar 
mass should satisfy m ^> H eq , where H eq is the value 
of the Hubble parameter at the time of radiation and 
matter equality. If we denote by t* the time at which the 
scalar mass equals the Hubble parameter m = H*, then 



772 2 = 



8vr 
3M] 



Pr- 



(5) 



We are assuming here that a large part of the scalar en- 
ergy density p^ was transformed into radiation, dominat- 
ing the matter content of the universe, and the rest is still 
in the form of an oscillating field. For t > t* the scalar 
field <f> oscillate around the minimum of the potential. 
The energy density average behaves as dust 



and also the radiation component evolves 



PR = Pr 



V a 



(6) 



(7) 



What we want to compare with observations is the dark 
matter mass per photon ratio ^ m = p^/n^. This quan- 
tity is constant for t > t* apart from changes in the num- 
ber of degrees of freedom of the species. We assume here 
an adiabatic expansion where S — gsT 3 a 3 is constant 
during the evolution, where gs is the entropic degrees 
of freedom, usually very similar to the total degrees of 
freedom g*. Because n 7 = 2((3)T 3 /n 2 , then 
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where the zero subscript indicates current values. Notice 
that we are assuming a change in the number of entropic 
degrees of freedom. Then the current dark matter per 
photon ratio is 
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The observational measure of this ratio is ^m.o = 2.2 x 
10 _28 M p using values from WMAP3 [9| , which for typical 
values ~ 100, gs.o — 3.9 gives the following constraint 
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4 x 10 
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(10) 



Using considerations from structure formation, we get an 
upper bound for the scalar field mass m/M„ = 10~ 52 or 
to > 10~ 23 eV. Obtaining the correct amplitude for scalar 
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perturbations requires m/M p ~ 10 -6 , then the condition 
(TIT)]) imposes that ~ 10~ 13 M p which is in contradic- 
tion with the initial statement for a value of <fi ~ M p at 
the end of inflation. In this case to get an observable vi- 
able model we need an incomplete reduction of the scalar 
field amplitude during reheating [l[ , reducing the energy 
density a factor of 10 26 . 

As I stressed at the end of section II, even in the stan- 
dard picture of reheating we can have a partial decay 
of the inflaton field. In this case, we find that at the 
beginning of the oscillations the field amplitude is al- 
ready ~ 10 _2 M P . So, the constraint mentioned in the 
last paragraph implies that during reheating the energy 
density must decays in 22 orders of magnitude. 



V. REHEATING AND THERMAL MASSES 

As we can see below, both features are connected; the 
presence of radiation during the perturbative reheating 
phase enable us to stop the process. Let us assume that 
during the process of preheating the inflaton decay prod- 
ucts scatter and thcrmalizc to form a thermal background 
Q. This thermalized particle species acquires a plasma 
mass m p (T) of the order of vT where v is the typical cou- 
pling governing the particle interaction. The presence of 
thermal masses imply that the inflaton zero mode cannot 
decay into light states if its mass m is smaller than about 
vT . If we expressed the inflaton zero-temperature decay 
width as = asm, at finite temperature it becomes 



IV. THE PROBLEM 

So, what we need, is to have an incomplete reheating 
phase after inflation, in agreement with the observational 
constraint: that of having at the beginning of the oscilla- 
tions (which is identified with dark matter) an scalar field 
amplitude 0* ~ 10~ 13 Af p . Clearly for a model where the 
amplitude of oscillations is very large, we can not neglect 
the parametric resonant effects during reheating Q. 

Using preheating, the authors of [l[ found that, for an 
interaction term g 2 x 2 4' 2 where g is the coupling constant, 
the decay is incomplete, stoping until the amplitude of 
the scalar field falls below m/g. Using the required ampli- 
tude of the scalar field at the end of inflation (see ea.(|10[)). 
0* ~ 10~ 7 m implies a coupling constant value of g ~ 10 7 , 
which is clearly incompatible with the model. 

However the restriction can not be applied at the end 
of preheating, because after the broad and narrow para- 
metric resonance phases Q , the reheating process follows 
through the standard mechanism. So, these conditions 
have to be applied after the entire reheating phase ends. 
In particular, the fact that the amplitude falls below m/g 
only means that resonant production stops. 

For a range g ~ 10 _1 — 10 -3 of the coupling con- 
stant, the amplitude of the field after preheating is 
$ ~ 10" 5 Aip-10" 3 M p respectively. These are the initial 
values for the standard reheating phase, so the observa- 
tional constraint implies that during particle decay the 
energy density has to decrease 16 orders of magnitude 
and not 26 as was settled in [l[. 

So, to get a viable observational model we have to con- 
sider the process of reheating once the resonant phase 
has finished. In this new context, the study of reheating 
is clearly different from the standard scenario; now the 
universe not only has the contribution of the coherent 
oscillations of the inflaton field but also the particles cre- 
ated during preheating. Therefore, our problem is to find 
a mechanism that not only allows to avoid that the field 
inflaton decays completely, but also take into account the 
presence of particles created in the previous phase. 



r^(T) = a^my 1 



2 T 2 



The system of equations to be solved is then 
P0 + (3H + r^)p0 = 0, 

p R + {ah p R - r^pj, = o. 



(ii) 



(12) 



When the plasma mass m p {T) = vT becomes comparable 
to the inflaton mass m, the temperature reaches the value 
T ~ m/2v, remaining constant for a while, indicating 
that particle creation stopped; = 0. At this time 
Pr stays constant and p<p decays as a~ 3 . In the absence 
of plasma masses, right after inflation the temperature 
increases until T max , and then decreases as a 3 / 8 until it 
reaches the reheating temperature T r h fioj 



100 
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T rh = 0.2 [ ) a 1 l 2 ,/m~M~ p . 



(13) 



In the case m < T r h, and using the values m ~ 10 8 GeV 
and g» ~ 100 (the effective number of relativistic degrees 
of freedom) we obtain 



> 3 x 10 



-10 



(14) 



To estimate the duration of this T constant phase, I use 
the fact that decays are not possible for T > m/2g. The 



scalar field energy density 



,-3 



until it becomes 



smaller than pr. So the condition is simply 

vf^y 4 



(15) 



where a* indicates the end of inflation and a/ the time 
when the <f> energy density equals pr. So, the universe 
evolves 



a f 



V 



1/3 



(16) 



which is equivalent to a number of e-foldings N e = 14 
after the standard particle decay process started, using 
the values V 1 ' 4 - 10 13 GeV and m ~ 10 8 MeV M. After 
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PR > P<t>i the universe becomes radiation dominated and 
the rest of the inflaton energy density p^ continue its 
evolution as a coherently oscillating field, behaving as 
dark matter. We can estimate the order of magnitude of 
the decrease in energy density using (fT6)) 

p^pi^j ^p^\ (17) 

which is in perfect agreement with the calculations of 
section IV. 



VI. SUMMARY 

In this paper, I have investigated the possibility to 
build a theoretical model which describe both, the infla- 
tionary phase and also the dark matter, using the inflaton 
field. The crucial ingredient is the period of reheating af- 
ter inflation, where almost all the energy stored in the 
inflaton field is converted into relativistic particles. In 
fact, as the authors of [1| have stressed, the main con- 



dition that must be satisfied is that not all the energy 
density of the inflaton decays into radiation. I found 
that, even in the standard (perturbative) reheating sce- 
nario, the process may end before all the inflaton energy 
density transformed into radiation. However, this obser- 
vation does not help to build the model by itself, because 
for a large initial amplitude of the inflaton field (with a 
typical value <pi ~ M p at the end of inflation), we can 
not neglect the parametric resonance effects Q. Taking 
into account the entire reheating phase, starting with the 
broad and narrow resonance phases, and after that, the 
standard perturbative one, I found that this can be ac- 
complished. The key element in this description is that 
inflaton decay products acquire plasma masses [|[ , which 
may cause inflaton decay to be kinematically forbidden, 
stopping the process. 
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